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Abstract 

We describe a high statistics quenched lattice QCD calculation of the 
moments of the polarized deep-inelastic structure functions gi and g2 of 
the proton and neutron. 



1 Introduction 

The polarized nucleon structure functions gi and g2 have attracted a lot of in- 
terest in the last years. Especially the measurements of gi [1] provoked a lively 
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discussion under the headline "proton spin crisis". The structure function g2, 
on the other hand, offers the first possibility of a direct measurement of a higher 
twist (=3) operator matrix element [2]. Experiments that perform measurements 
of g2 are being undertaken at DESY (HERA) and SLAC. On the theoretical side 
it is widely believed that quantum chromodynamics (QCD) is the correct theory 
of the strong interactions. To put this belief on a solid experimental footing one 
has to compare the experimental results with the predictions from QCD. It is 
therefore a major challenge for theoretical particle physics to find answers to the 
question "What does QCD tell us about the structure functions in deep-inelastic 
lepton nucleon scattering?" 

Asymptotic freedom allows us to calculate the violations of Bjorken scaling in 
the (moments of the) structure functions using perturbation theory and Wilson's 
operator product expansion (OPE). However, the computation of the structure 
functions themselves (for a fixed value of the momentum transfer Q^) requires a 
nonperturbative method. This talk describes an attempt [3] to perform such a 
calculation applying lattice gauge theory and Monte Carlo simulations. (For ear- 
lier work see ref. [4].) We have studied polarized as well as unpolarized structure 
functions, but in accordance with the subject of the workshop we shall restrict 
ourselves to the polarized case in this talk. 

From the OPE we have for the leading twist contributions in the chiral limit: 



f=u,d 

Jo Zn + I f_, L 



for even n and n > (n > 2) for gi (5^2 )■ The reduced matrix elements a[{\ d[l^ 
for the flavors f = u,d are defined in terms of forward nucleon matrix elements 
of the operators 

^S--.. = (^'^ ^7.75 D,,--- D,^i; - traces (2) 
with -0 = u{d) for / = u{d) according to (see e.g. ref. [2]) 

(P^ ^\^{^Z■■■^^n}\P^ ^) = ^^[^-^'''1 ■■■P^^n + traces] , (3) 

5(f) d''^'' 
(P^ ^P[a{nr]-nn}\P' ^) = ^^^[(^^^''i ~ Sf,,p^)pf,, ■ ■ -Pf,^ + traces]. (4) 

Here {• • •} and [• • •] denote symmetrization and antisymmetrization, respectively, 
and fi is the renormalization scale at which the operators are defined. In the fiavor 
singlet case purely gluonic operators will contribute also. 
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The Wilson coefficients ell{ii^ /Q"^ ,g{ii)), j = 1,2, are calculated in pertur- 
bation theory. For electroproduction they have the form 

e^%yQ^g{^i)) = Q^f^Hl + g{f^)%n{f^yQ',gm , (5) 

where Q^-^^ are the quark charges. The nucleon states \p, s) with momentum p 
and covariant spin vector s are normalized according to 

(p, s\p', s') = {2Trf ■ 2E^d{p- p')Ss,s' (6) 

where = —m% {rriN = nucleon mass). Note that C'll^^^...^^^} has twist 2, whereas 

^[a{lii]-n„} '^^^ operator of twist 3. Hence, in the moments of ^'2, clU^ represents 
the so-called Wandzura-Wilczek contribution [5]. 

The basic task is now the computation of the matrix elements (3), (4) in a 
polarized nucleon. Depending on the choice of the indices cr, /ii, . . . , //„ it will be 
necessary to give the nucleon a nonvanishing three-momentum. 



2 Lattice Calculation 

How are these matrix elements calculated by lattice Monte Carlo simulations? 

Analytic continuation to imaginary time, i.e. putting xo = —1x4, leads us from 
Minkowski space to Euclidean space and allows us to interpret the functional inte- 
gral of a quantum field theory as the ensemble average of a statistical mechanical 
system in four dimensions. We can then take advantage of the experience gained 
in statistical mechanics, in particular we may discretize the Euclidean space and 
use Monte Carlo simulations for the calculation of the expectation values we are 
interested in. 

Roughly speaking, the continuation to imaginary time means that the unitary 
time evolution operator e~'^* is replaced by the self-adjoint operator e~^*. Hence 
the exponential decay of Euclidean correlation functions, the so-called Schwinger 
functions, contains information about the energy spectrum of the theory, and the 
amplitudes multiplying the exponentials are related to matrix elements. 

To be more specific, consider the two-point function 

{0\B{t2)B{ti)\0) = (0|e'^*=5(0)e-'^*=e'^*iB(0)e-'^*i|0) 
= (0|S(0)e-'^(*^-*i)E(0)|0) 

with suitable interpolating fields B, B for the nucleon. Analytic continuation 
yields the correlation function 

(0|Be-^*5|0) = {B{t)B{0)) . (8) 
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As t becomes large, it is dominated by the lowest-lying state coupling to B, i.e. 
the nucleon with energy E^: 



{B{t)B{0)) = {0\B\N)e-^^\N\B\0) + ■■■ (9) 

In order to suppress the contributions of higher excitations indicated by the dots 
in (9) we use the freedom in the choice of the interpolating field B by applying, 
e.g., 'smearing' techniques [3]. 

For three-point functions we have similarly 

{0\B{t)O{T)B{0)\0) = (0|B(0)e-'^(*-")e'(0)e-'^^B(0)|0) (10) 

In Euclidean space we obtain for t > r the correlation function 

(0|Be-^(*-")Oe-^^5|0) = {B{t)O{T)B{0)) (11) 

and in the limit of large time differences t — r, r we have 



{B{t)O{T)B{0)) = {0\B\N)e-'^''^'-^\N\O\N)e-'^''^{N\B\0) + 
= {0\B\N)e-^^^{N\B\0){N\O\N) + • • • 



(12) 



Hence the leading contribution is independent of r and we can determine the 
desired nucleon matrix elements from ratios of three-point functions to two-point 
functions: 

<^|||^.„)......0<<.<<.. (.3, 

In order to regularize the ultraviolet divergences we introduce a regular hy- 
percubic lattice with lattice spacing a. Our lattice has the size 16^ x 32, with 
the larger extent corresponding to (Euclidean) time. The quark fields are chosen 
periodic in space and antiperiodic in time, whereas the gauge fields are periodic 
in all directions. We take the standard Wilson action for the gauge fields and 
use Wilson fermions to put the quarks on the lattice. The latter choice implies 
an explicit breaking of cliiral symmetry by lattice artifacts and the bare quark 
mass is replaced by the so-called hopping parameter k. The quark mass is then 
approximately given by 

m„a ^ 0.56 (- - —] , (14) 

where the critical hopping parameter Kc corresponds to the chiral limit (vanishing 
pion mass). Our bare gauge coupling is g = I and we have data for three values 
of the hopping parameter, k = 0.155, 0.153 and 0.1515, corresponding to quark 
masses of roughly 70, 130 and 190 MeV, respectively. Hence we can try to 
extrapolate our results to physical quark masses. For the accuracy we are aiming 
at this is the same as extrapolating to the chiral limit. 
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Our calculations are performed in the so-called "quenched" or "valence" ap- 
proximation, i.e. we neglect all internal quark loops, as this drastically reduces the 
necessary computer time. Correspondingly, we disregard in the operator inser- 
tion the so-called "disconnected" contribution where the operator is not directly 
coupled to the valence quarks in the nucleon. We are, however, trying to calcu- 
late this contribution thus making a first step towards incorporating sea quark 
effects, but we have no results yet. 

3 Operators and Matrix Elements 

In Euclidean space we deal with operators of the form 

^Im-nn = "^T^Ts^Mi • • • ^Mn^ - traces . (15) 

On the lattice, Df, is of course the lattice covariant derivative. More explicitly, 
we study the three operators 

O2 = -07275-0 forao , (16) 
(^{214} = ^ + Ol, + Ol,, + ol,, + ol,, + ol,,) for a, , (17) 

^f2{l]4} = 2 (^214 + ^241 ~ 1^124 ~ 1^142 ~ 1^412 ~ 1^42l) ^^1 d,. (18) 

Two criteria guided us towards the above choice of operators. First, we do not 
want to be forced to have more than one component of p nonzero. For this would 
lead to higher energies of the nucleon states to be considered, hence to a faster 
decrease of the propagator and ultimately to a worse signal-to-noise ratio. 

Secondly, we have to avoid mixing with lower-dimensional operators. These 
would come with coefficients diverging like a power of a^^ . So the necessary 
subtractions would be large and furthermore uncalculable in perturbation theory 
[6]. Note that the mixing pattern is much more complicated on the lattice than 
it is in the continuum, because the constraints imposed by symmetry are less 
stringent. Instead of covariance with respect to the orthogonal group 0(4) we only 
have at our disposal the finite subgroup H(4) (the so-called hypercubic group) 
which leaves the lattice invariant. Multiplets that are irreducible with respect to 
0(4) may become reducible under H(4), and the possibilities for mixing increase. 

Corresponding to the above operators we choose the polarization in the 2- 
direction and consider the two momenta p = , (27r/16, 0. 0) • 1/a. Furthermore, 
we have to choose a fixed value for t in eq.(13). This should be as large as possible 
in order to allow for a long "plateau", i.e. an interval in r where the ratio (13) 
is independent of r and hence enables a reliable determination of the desired 
matrix element. We have taken t = 13. The ratios pertaining to the reduced 
matrix elements ao, a,, and d, are shown in Fig.l. 
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Figure 1: The ratio R for u and d 
quark insertions for k = 0.153. (a) 
-ii?ao, (b) Ra2 and (c) 2- T IS 
given in lattice units. The horizon- 
tal lines indicate the fitted value of 
R on the plateau together with its 
error and the fit range. 



10 12 14 16 18 20 



(c) 



For the three operators (16)-(18) the ratios R of three-point functions over 
two-point functions (in lattice units) can be written as 



■ao 



R. 



•02 



R 



•d2 



1 1 



Zao 2/« 2 



Qq for p = 6 . 



1 1 rriNPi 



Za2 2/« 



6 

J- 1 rriNPi 

2k 3 



0-2 , 



do 



(19) 
(20) 
(21) 
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The factor 2k results from the different normalization of the lattice quark fields 
as compared to the continuum conventions. The Z's {Z = Z{aii)) are the renor- 
malization factors needed to construct operators renormalized at a scale jj, from 
the bare lattice operators. At present, we take the Z's defined according to the 
momentum subtraction scheme from a one-loop perturbative calculation in the 
chiral limit. (A nonperturbative calculation seems however possible [7] and is 
being done.) If necessary, we can convert the numbers to the MS scheme, in 
which the Wilson coeflBcients are usually calculated. In the following we shall 
quote our results extrapolated to the chiral limit for 

^ = Q^ = ^ 5GeV^ , (22) 

which eliminates the logarithms in the Wilson coefficients and renormalization 
constants. 

Let us note in passing that the operator (18) used for the determination of 
^2 could mix on the lattice with some other operator of the same dimension. 
However, the matrix element of the latter operator turns out to be small and 
the mixing coefficient vanishes in one-loop perturbation theory. Hence we have 
neglected the mixing. 



4 Discussion of Results 

In Fig. 2 we show our results and their extrapolation to the chiral limit. Instead 
of Qq we plot A?i and Ad defined by 

= 2Am , a-i^ = 2 Ad (23) 

and interpreted as the fraction of the nucleon spin that is carried by the u and d 
quarks, respectively. For the proton we find in the mom scheme 

Au = 0.83(7) , Ad = -0.24(2) . (24) 

Since these numbers were obtained in the quenched approximation, they should 
be considered as an estimate of the valence quark contribution only. A quantity 
which is expected to be less sensitive to sea quark effects is the difference Au — Ad 
for which we get in the proton: 

Am- Ad = 1.07(9). (25) 

This is to be compared with the experimental value of the axial vector coupling 
constant of the nucleon, which is 1.26. 

Instead of looking for observables which are dominated by the valence quarks 
one can try to incorporate at least part of the sea quark effects by taking into 
account the "disconnected" contribution to the matrix element. Our calculations 
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Figure 2: Aq, a2, and c?2 for the 
proton as a function of together 
with a linear fit to the data. The 
solid symbols indicate the extrap- 
olation to the chiral limit. In (a) 
we compare our numbers with the 
phenomenological values of Ref.[9], 
which are marked by asterisks. (The 
value for Au is hidden behind the 
lattice number.) 



of this quantity are not yet finished, but for Aq we can use the results of two 
other groups [8]. Renormalizing these results according to our conventions we 
obtain in total 

Au Ad As AS 

Lattice QCD 0.69(9) -0.38(5) -0.13(4) 0.18(11) 
EK 0.83(3) -0.43(3) -0.10(3) 0.31(07) 

For comparison we show in the row labeled EK the results of the recent phe- 
nomenological analysis by Ellis and Karliner [9] using = lOGeV^. 
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Rewriting our results in terms of the lowest moment of gi we get (cf. eq.(l)) 



/' dx 01 fx O^) - I 0-166(16) proton, 

Jo "^^^U^'^ ) - \ -0.008(09) neutron, ^^^^ 

with the "disconnected" contributions left out. This omission should be less 
harmful in the difference of proton and neutron structure functions. Indeed, our 
result 

/' dx{ff,{x, Q') - 9^{x, Q')) = 0.174(25) (27) 

compares rather favorably with the experimental value of 0.161 ±0.007±0. 015 [9]. 

In the higher moments of gi sea quark effects are expected to be suppressed. 
From our calculation of a2 we obtain in the chiral limit 

f'dxx'oAx 0')-l 0-0150(32) proton, 
Jo <^^^ 9i[x,LJ ) - <^ -0.0012(20) neutron, 

consistent with experiment. 

Let us finally come to the structure function g2. Extrapolating our results for 
a2 and ^2 to the chiral limit we get 



^0 



'•1 2 1 _/ -0.0161(16) -0.0100(22) = -0.0261(38) proton, 

/o f2(2;,y ) - I .QQQig^Qgj +0.0009(13) = -0.0004(22) neutron. 

(29) 

Here the first number comes from ^2, while the second number comes from 02. 
Hence the twist-three operator provides a significant contribution, at least in the 
case of the proton. This result disagrees with sum rule calculations [10], which 
suggest that for the proton ^2 is very small. 



5 Conclusion 

We have presented a calculation of moments of the polarized nucleon structure 
functions from first principles. The results show that such a calculation is feasible, 
at least for the lower moments. We had to make some approximations though, 
e.g. we have neglected sea quarks and gluonic operators. Furthermore, there 
are systematic uncertainties due to the extrapolation to physical quark masses, 
cutoff effects etc. However, all these shortcomings can (and will) be improved. 
But even at the present stage the agreement with the experimental data, as far as 
they exist, is rather satisfactory. More precise experiments measuring g2 should 
reveal if the twist-three contribution indeed plays such an important role as our 
results suggest. 
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